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Abstract 

The hydrodynamic conservation equations and constitutive relations for a bi¬ 
nary granular mixture composed of smooth, nearly elastic, hard spheres with 
non-equipartition energies and different mean velocities are derived. This re¬ 
search is aimed to build three-dimensional kinetic theory to characterize the 
behaviors of two species of particles suffering different forces. The standard 
Enskog method is employed assuming a Maxwell velocity distribution for each 
species of particles. The collision components of the stress tensor and the other 
parameters are calculated from the zeroth- and first-order approximation. Our 
results demonstrate that three factors, namely the ratios between two granular 
masses, temperatures and mean velocities all play important roles in the stress- 
strain relation of the binary mixture. The collision frequency and the solid 
viscosity escalate with increasing of two granular temperatures and are maxi¬ 
mized when both of two granular temperatures reach maximums. The hrst part 
of the energy source varies greatly with the mean velocities of spheres of two 
species, and further, it reaches maximum at the maximum of relative velocity 
between two mean velocities of spheres of two species. 
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1. Introduction 


Granular mixtures [DEIIS], including e.g., landslides, avalanches and phar¬ 
maceutical powders, are pervasive in nature and industry. Granular mixtures 
and ordinary molecular gas mixtures act differently. For a granular mixture 
5 particularly engaged in rapid flow, nearly instantaneous and inelastic collisions 
happen between particles, leading to dissipations. So it is necessary to inject 
energy continuously to keep a steady state for a granular mixture. Two ways 
of energy compensation are usually involved: boundary driving, as vibration 
[aa] and shear [5] , and bulk driving, as air-fluidization iE] and magnetic field 
10 [S]- Similar to a single species granular gas, granular mixtures are normally in 

non-equilibrium states. And it has been demonstrated that the components in a 
granular mixture do not share the same granular temperature [a uni EH [la EH- 
The dense gas kinetic theory of Chapman m has been employed to quan¬ 
tify rapid flow of granular mixture systems for decades E3 El EH EH], in which 
15 granular mixtures are assumed to be smooth, nearly elastic and spherical grains. 
In literature, two approaches can be classified, as elaborated by Galvin et al. 
m- the first one is derived via systematic expansion, such as Chapman-Enskog 
method, and the second one is on the base of a hypothesized velocity distribu¬ 
tion. According to the factors of the velocity distribution (Maxwellian or non- 
20 Maxwellian velocity distribution), equipartition or non-equipartition (equal or 
unequal granular temperature), radial distribution(standard Enskog or revised 
Enskog theory [ID]), system dimension (two or three dimensions), the second 
approach can be further divided. Jenkins and Mancini [20] firstly derived consti¬ 
tutive relations and balance equations by employing the Enskog equations with 
25 the assumption that the velocities of two species of particles are Maxwellian 
distributed, and the temperatures of two species were supposed to differ by in¬ 
finitely small quantity. Based on the revised Enskog theory [211 EH EH] , more 
exact theories with equipartition assumption were further developed, concerning 
the non-uniform, local equilibrium sate. The kinetic theory of non-equipartition, 
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30 binary granular mixtures with different sizes, masses and diameters was also 
obtained based on the Enskog equations by Huilin et al. [SUES], where two 
granular temperatures were assumed, respectively, for different sized particles. 
However, both Rahaman et al. and Iddir & Arastoopour [55] pointed out 
that Huilin’s collision rate between particles i and j, Nij, is not symmetric, 
35 that is, Nij ^ Nji. Rahaman et al. [25] improved the integration processes by 
supposing the angle between the relative and combined velocities of two species 
of particles i and j is in the range of [0, 27r], implying both of these velocities 
are two-dimensional vectors. That treatment contradicts the three-dimensional 
derivation in their work. Moreover, all the above derivations supposed that the 
40 mean velocities of two species are identical, which is obviously not the case in 
gas-fluidized systems |571[li[21[3D]. The first approach developed by Garzo et 
al. |3TJ |32| was dedicated to solving perturbatively the Boltzmann equations to 
capture a broader range of restitution. However, Galvin et al. m pointed out 
that the application of the second approach could be extended to moderately 
45 dense system, whereas the first one is limited to dilute granular gases. Recently 
Garzo et al. [33l |34] extended the first approach to the (moderately) dense 
flow. Both of these approaches have been used to model the stress of particles 
in multiphase computational fluid dynamics (CFD) [HI [5?11^ ITfl IM| . 

It is worth noting that not only the binary granular mixtures but the air- 
50 fluidized granular systems with identical spheres are also in non-equilibrium 
states pni Ell EH EH- III gas-solid circulating fluidized beds, dense clusters 
El HD] with size of 10-100 times the particle diameter are suspended in the di¬ 
lute broth of gas-solid mixture, whereas the particles accumulated in the cluster 
and those dispersed in the dilute broth do not share the granular temperature 
55 and mean solid velocity [10] ■ This is reasonable because the particles in the 
dense cluster and dilute both suffer different drag forces [271130] which make 
the particles away from equilibrium states (share the same granular tempera¬ 
ture). To illustrate this inhomogeneous structure, it is convenient to elucidate 
the above situation by using two velocity distributions with two temperature 
60 equations. Thus, the kinetic theory of double granular temperatures and mean 
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velocities makes sense not only for the binary granular mixtures but also for 
CFD simulation of fluidization. 

In this paper we dedicate to deriving the constitutive relations and bal¬ 
ance equations of two kinds of particles with unequal mean velocities and non¬ 
es equipartition energies using the standard Enskog theory. Sec. provides a 
binary collision frequency derived through their dependence on the velocity 
distribution function of a mixture of inelastic, hard spheres. In Sec. the 
hydrodynamic descriptions are given based on the Boltzmann equation. In Sec. 
1 ^ the constitutive equations are identified by macroscopic hydrodynamic vari- 
70 ables. In Sec. we discuss our results and compare with previous works. In 
Sec. the main conclusions are summarized. 


2. Binary collision frequency 


For a binary granular mixture composed of hard, smooth, inelastic spheres 
of species a and /3 with mass and diameter di, i = a,/3, due to suffering 
inhomogeneous external energy input, such as the different drag forces exerted 
on fluidized particles, two species may have their own temperatures and mean 
velocities, respectively. Here, we consider the most general case for binary 
collisions. Let i and j represent either species a or f3. For an inelastic collision 
between two particles in species i and j, with velocity c^ and Cj, the relationship 
between pre- and post- collision relative velocities c^- = c^ — Cj and cF = c' — c' 
yields 


k • c^j (k • Cjj) 


( 1 ) 


where k = Vij/rij is the unit vector directing from the center of particle with 
velocity c^ to that of particle with velocity Cj upon contact, specifying the 
geometry of the impact. Cij is the restitution coefficient between species i and 
j. Angular velocities or rotations are not included in this paper. The velocity 
of the center of mass G is defined as: 


G = 


rrii Ci + rrij Cj 
mo 


( 2 ) 
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where ttiq = rrii + rrij. At the moment of collision, the distance between two 
centers of particles is dij = {di + dj) 12. The number of binary collisions between 
species i and j at position r per unit time per unit volume has the form: 

Nij = / fj;f {ci, r, Cj,r + dyk) (cy • k) dj^dMcidCj (3) 

J Cij ■k>0 

( 2 ) 

where ' is the pair distribution function. Following Chapman and Cowl¬ 
ing El- the assumption of chaos allows us to write the correlation probability 
function as a product of two single velocity distributions: 

(ci, r, Cj,Y + dyk) = x(r -h ^d*jk)/, (c,, r) fj (c^-, r dyk) (4) 


where the factor y called Enskog factor equals to unity for a rare gas [14], and 
reads m- 


X(r -h -dyk) = gij 


9ii ^j') “f ^i9jj ) 
2dy 


( 5 ) 


in which 


9ii 5 ) 


1 


1 — (Ci + £j)/£ 


max 


3di £k 

2 dfc 

k=i,2 


( 6 ) 


And £i and £j are solid volume fraction for species i and j, respectively, £max 
is the single-phase maximum packing, gij is the radial distribution function 
75 between spheres of species i and j. 

And we assume the velocities of both species of particles follow the Maxwellian 
distribution: 


f^ (ci, r) = n, 


f 


\2TT9i 


3/2 


exp 


rrii {ci - v^y 


29,. 


( 7 ) 


where rii is the particle number density, is the mean velocity of species i, 
and 9i is the granular temperature defined as an ensemble average 9i = ^rrii < 

(Ci - Vi)2 >. 

For the binary mixture, the mean mass center velocity is defined by 


niiTiiVi + njiTijVj 
niTTii + njrrij 


(8) 
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The peculiar velocity is the relative velocity between the particle velocity and 
the mean velocity 


The diffusion velocity reads 

Uj = Vj - V 


(9) 


( 10 ) 


which denotes the mean velocity of species i relative to the mean mass center 
80 velocity. 

Then the joint pair distribution function becomes 


r( 2 ) f \ 1 

fij (Cj) Ti, Cj, Tj) — 


/ \ 3/2 

f miTrii \ 

TOi (Ci - Vif 

1- 

CM 

> 

1 

[ 0,9, ) 

1 

1 

to 

1 

20, 


( 11 ) 


Using this pair distribution, the collision frequency becomes 
1 / \ 3/2 

AT 1 72 ! ITT-imj ' 




Cij exp 




20 ,, 


26, 


( 12 ) 


dcidcj 


We can express the particles velocities in terms of G, Cj, as in literature 
El [^. by expanding it in a Taylor series: 

3/2 


1 

Ni, =—^d. 




miiTij 


Cij exp [-AG^ - Dc 




(13) 


X [1 - 2B(G • c„) + 2B‘^{G ■ c,^)^ + ...]dGdc 


where 


^ ^ m,6j + rujOj ^ ^ 

26,0j ’ 


mimj{6i - 9j) 
2731^6 id j 


D = 


mimj{mi6i + rujOj) 
2mQ6i9j 
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cos(a:) = cos(^i)cos(^2) + sin((Pi)sin(^2)<^os(0i “^2) 


Figure 1: The sketch of a binary collision. 


The integration of the cross term G-Cy needs to be stressed here. In Huilin’s 
model [21] , the mass center velocity G and relative velocity c^- in the cross term 
are treated to be scalar. In Rahaman’s model [25] . G and c^- are assumed in a 
two-dimensional plane, so the angle between G and Cij are in the range of [0,27r]. 
However, as illustrated in Fig. 0. in three-dimensional spherical coordinates, 
the angle between two vectors G{ri,(pi,9i), c^-(r 2 , (^ 2 , ^* 2 ) is 


a = arccos(cos ipi cos ip 2 + sin(/Ji sin ip 2 cos(0i — ^ 2 )) 


(14) 


Integrate Eq. (13) using Eq. (14), we obtain 


Nij = -g.jdtjn.Tij 


/ rriimj ^ 

N ^ f 

V 

/ A^ND'^ 


456 


AD A^D^ 


2B^ 3B 

1 + -TTV + 


(15) 


Due to the symmetry of collisions, the number of collisions of species i is 
expected to equal the number of collisions of species j, i.e., Nij = Nji, which im¬ 
plies that B in the expression of Nij should be raised to power of even numbers. 
85 Our results satisfy this criterion. 
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3. Conservation equations 


The velocity distribution function / (c, r, t) of each granular species of a 
binary mixture subjected to an external force satisfies the Boltzmann (integro- 
differential) equation [T3] 



(16) 


where ( ) is the change rate of distribution function / (c, r, t) due to par- 

tide collisions, F represents the external force such as gravity, buoyancy and 
gas-solid drag force. 

For species i, let ijji be any function of particle velocity c^. Multiply both 
sides of Boltzmann equation by tpidci and integrate them throughout the velocity- 
space, we could obtain the equation of change rate of the particle property, the 
transport equation for the quantity '0i: 



drii <ipi> 

dt 


-|-V • I < *-2^02 > T ^ ( ^c,ij 


(17) 


where Pc,ij is the collisional part of the stress tensor, and the collisional 

source term for particles in species i during the collision with particles in species 
j. Here, we use the expressions of Pc,ij and given by Jenkins and Mancini 
EOl in 1987: 



(18) 



k Cij >0 


(0' - 02) (k • Cij) 



(19) 


Setting ipi = rrii, we get the transport equation of mass. 



( 20 ) 
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Setting ipi = rriiCi, we get the transport equation of momentum. 


d{nimiVi) 


dt 


- V • 


+ V • (n^miViVi) = 




+ nimjFi+ ^ 


( 21 ) 


Among it, is the kinetic part of the stress tensor for species i. 

Pfc.i = J m^CiCifidci (22) 

Setting ipi as ^micf and subtracting it by the product of the momentum bal¬ 
ance equation (Eq. ( |21[ )) and hydrodynamic velocity v^, we obtain the equation 
of granular temperature, 


■ d 

— (nj 6 »i) -t V • {riiVie,) 

= 




II 


-t- TliTTli < F • Cj > -j- En., 


: Vvi - y{qk,i + (lc,t) 
1 






(23) 


where j represents J2j=a p ^c,ij ( 5 ^ 1 ^^), and is the kinetic part of en¬ 
ergy flux which is given by 


q/c.i = / ^m.CiCJidc, 


(24) 


The above hydrodynamic equations (Eq. (|20[), Eq. (21), Eq. (23)) are based 


on the continuity hypothesis for the granular system. So the established range 
of our results are limited to the Navier-Stokes level. 


4. Constitutional relations 

In this section, we aim at deriving several terms, such as Pk,i and q^. j. In 


95 these derivations, the pair distribution function is employed as Eq. (11 1 . 


To begin with, we calculate the stress tensor P which is caused by two 
mechanisms, each contributing one part to P [H]. The first is the kinetic part 
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of the stress tensor, denoted by Vk^i as defined in Eq. (22). It results from 


the motion of all particles with mean velocity without any effect of particle- 
100 particle collisions. 

Calculated from the first-order approximation to the distribution fi, Pk,i 
becomes 


Pfc.i = n^dii - M + + e^) ) V^u,- 

gu \ 15 


(25) 


where I is the unit tensor, is the restitution coefficient between particles in 
species i, and gi^du is the viscosity for dilute suspensions [5], expressed as 




hrm-J-KOi 


(26) 


167rdf 

The second part of the stress tensor P describes the momentum transfer 
caused by collisions. It is denoted by Pc,y(miCi). We obtain it by setting 


V'i = m^Ci in Eq. (18). 


— o (k ■ ) k 

^ ./kCn>0 


X ( Cj, r - k, c^-, r -h ^dijk ) dkdcidcj 


(27) 


We still adopt the expansion of (ci,r — Id^jk, Cj, r-|- ^dyk) in the paper 
[23], then, 


Pc,ij(^iCi) ^ ^i)(^b • k'^kfifjdkdc^dcj 

j ^ dci,k>0 

- yffijdfj f TOj(C' - Ci){c,j ■ k)kf JjVln^dkdcidCj] 
4 dkcij>0 Jj 

(28) 

where Pl^j ^c,ij indicate above two integrations, separatively. We assume 
the diffusion velocity and u^- approaches zero, so a factor ej:p(——^u^) 
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becomes 1. Using the integration technology msa, we get 


I 1 + Cii rriim^ 


P = 

c,n 


-ij ,3 


48 mo 


I mi TTlj 


3/2 




, 5^2 ass'* 

1 + TTTT^ + 


2AD SA'^D^ 


(29) 


and 


= 

3847r2 


— BRoS 




j 

9i9j 


‘iji 


mi 

TTlj 

V 0 , 

V 0 U 


mi 

rrij 

mi's /9 i 

mjS/9j 

[ 0 . 

* e, 



+ 3526 

[0^ 

e, 

9^ 

J 


^ _B 

3mo mo 


'mi 

rrij 

m j V 9i 

<1 

+ 25i5 

Uj 

Uj 

S/9, 

V0,1 


9 , 

9! 

J 


9, 

9, 

J 


mfm| j 25 
m? I 3 


555 


-1 

c 

c 

<1 

h 

<1 

-^545 

Ui 


m j V 9i 

miS/9j 

9, 0/ 9] _ 

9, 

% 

91 

to 


miirij 


Ui Uj miS/9i mjS/9j 


mi mi S/9i S/9i 

1 

mo 

A 

1 

1 

_ 1 

+ 3 

9, 9j 9} ^ 9] 

J 


-1 


1 m,mi „ \ ~ ■, / 25 m,- „ 1 m,mi „ , ^ 

’■«.)5|V.v,l+(-^fl»-^^ft)5|V,v.l 


6j mo 


(30) 


where we denote a tensor S as the function of two vectors v and w, 


" 4 o 2 

5[v,w] = -[vw]'* + -[v • w]I 

t) O 


(31) 


vw means the dyadic of two vectors. For any tensor X, the superscript means 


i^ = i(x + x)-i 


(X:i) 


(32) 


where the notation X is conjugate of X. And the following substitutions are 


105 used in Eq. (301. 


4571^/2 _ 471^/2 ( 352 \ 1571^/2 _ 671^/2 ( 5 B ^\ 

° “ AV2£)4’ 1 “ A 3/253 AD y ’ ^ “ AV^D^ ’ ^ “ A 5/253 AD y ’ 
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127r5/2 / 4^2 \ _ 187r^/2 ^ 5^5/2 ^ g^5/2 / ^2 

4 - ^5/2£,4 +ADy’ 6 - ^5/2£,3’ ^ “ A5/2£,3 + ]^ 

Consequently, we acquire the total particle pressure Pj from the normal term 


of the sum of Pfe^i (Eq. (25)) and Pc,ij (Eq. (29), Eq. (30)) 


rriimi 1 , \ I "i-i" 

p, =n,e, + + l)dyg,,n,n, ( — 


Too 


niim-i 

w 


•ill 


1 


5^2 3554 

1 + + 


(33) 


A3I2JJ5/2 y 2AD 8A^D^ 

In addition, the coefficient of viscosity is defined from the the shear term 
of the sum of P^^i (Eq. (25)) and Pc,ij (Eq- (29), Eq. (30)), thus 


l^i dil f 27r o , , "N f 27r q , . 27r •n^TTi/'i q , , 

+ ei)j 1^1 + —(firi,gu{l + e,) + + ^ij) 


E 


1 / -| \ >4 / \ I 11 ^j 


QGOtt^ mo 


-(Cij + l)dijgijri^n 




/TOim, 1 2B m, 

Bi - — -pr^6 


mo 9i 3 9i 


(34) 


The momentum source term ipi involved in Eq. (21) is therefore 


Ny(miCi) 


= dy- [ {miC[ - rriiC,) (k • Cij) f[f ( c^, r - k, , r + E^jk ) dkdc.dcj 

dk-Cij>0 \ ^ ^ / 

= E 

3 

(35) 


where 


777 ■ TTl ■ / 

=-(1 + Cij) / (cy • kfkfifjdkdcidcj 

^0 7ci.j-k>0 


(36) 


^(1 + eij) f (cij •k)^k/i/jk-Vln^dkdcidcj (37) 

■'410 ^ dci^-k>0 Jj 
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110 whence, to the same approximation and integration, 


a _ (1 + Sy) mi'm-j ^2 „ „ „ ( ’niirrij 


12 


mo e,e 


3/2 




~D 


nii 


- B ( ^U,; - 


(38) 




Lr. 

V9A 


1 V rrij 


- Ui 


2 ^ /1 I \ j3 

= ~ ^ijjdijgijnirij 

+ [4(V • V,)/ + 8VvJ] • u, + 10 f 

\mi rrij 

20 , 


-u, • [4(V-v,)/ + 8Vv*] 


rrij ^ rrii^ 

^Vvj- • Uj- - ^Vv, • u, 

Uj Ui 


+ 8 {Vvj ■ u^- - V®Vi • Ui) + y [(^ • Vj)u^- - (V • Vj)ui] 


(39) 


As defined in Eq. (24), q^, ^ denotes the kinetic energy flux due to the 


particle transport without collisions. Using the zeroth-order approximation, 
qj, i is found to be equal to zero, since the integration is on the odd power of 
the velocity C^. 

Extending q,r. ^ to the first-order approximation, we obtain 


( 1 ) ^i,dil 

i = - 

gti 


1 + + ei)^ 




(40) 


115 where Ki^du is the conductivity for dilute particle phase, expressed as Hi^du = 
^gi,diiCv, and Cy = 2 ^/ 7 - 

The collisional contribution to the energy flux is determined by Pij (^m^C^), 
where 
'1 


7;m.Ct 


y ( 2' 

= “X 4 “ ^i) r - ^k, Cj, r -k ) dkdcidcj 


= - • 


(41) 


and 


rriiC, - 


1 




Cf) ■ c,j)kfifjdkdc,dCj (42) 
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2 ^ f 

4yk.c.,>oV2 


^m^Cf - (k • Cij) kfifjk ■ V In ^dkdcidcj 

(43) 


Among them, — Cf) is the kinetic energy change for species i in one 

collision. As noted above, the collisions between the particles are inelastic, so 
the kinetic energy change during a collision is 


1 


1 


C fZ ^ 

i = — 


^(1 + ■ k)^ - 2^(1 + ejj)(cy • k)(Ci • k) 


TOn 


Too 


(44) 


Substituting Eq. (44) into Eq. (42) and Eq. (431, we obtain 


qLi = (1 + eij)nin_ 


miTTij 




(u, - Uj) 


(45) 


2 TTliTTlj , ^ 


(^Ro + \rg 
mo \ 5 3 


V9i V0,- 


el 




- “^BR \ j _ 2^^ (mfmj V9j _ m^m^ V0, 


3" M 0? 0? 3 M ml 91 ml 9? 


/V9i V6»* 

^BReiY-ir 

Uj Ui 


2mi 


" (e.,-1) 


4 „„ /mimj V9i mimj 


3 V mo 9f 


+ 


I mjV9j miW9i i j, / "h"^3 

A-s 1 —- 7 ?)— I + -^4 


m^mj V9j m^m^ V0i 


0^3 


/V0, 

V0A 

'v % 

J 

(rrij ^ 

2B , 

-^Ri 

-q“-^ 

K^mo 

3 


''0 ''3 

2B „ 


9j V mo 


'jJ 


mo 91 


(46) 


Finally, we substitute ipt with energy l^miCf, then collision part of energy 


14 



























source term can be written as 


N.,( 


= 

/ 


-'k-Cij>0 

= 7^ - - + 7^ - - 

'C^lJ ' /C,2J 




7 ^ ^ 2 ^ ) (^^dcidcj 


(47) 


where 


7c,ij — d^jgij j ( ^rriiC^j (k. ■ Cij) fifjdk.dcidcj (48) 

^ k 'Ci 7 ^ 0 \ / 




^rUiCf - im*C-^ (k • ) /i/j V In y d^dcidcj (49) 


2 Jkcij>0 

Following the above integration method, ^c,ij includes 
1 1 M , n ,2 rrurrij fmimj 

— g ^3/2£)3 


mo 


6i9j 


3/2 



- I'l - 

2 mimj j 

^-Ui • Ui Ui • Uj V 

_mo 

-•-j 

3 mo ' 

V J\ 


2B 

'si 


^ j-^ I TTli no'i 

+ n - 1 3 - -^u, • Ui - —Ui ■ Uj 


1 + 


~AD ) 


(50) 


and 


2 TT miTTlj ^ 2 

^cAj ,, 2 

’ 4 TOq 


(cjj l) d^jgijTiirij 


m, m-i 


V6I, 


- + V.v. + V.v, + ^.u 


9. 


S/9,, 


+ - 

9^ 


f 1 

X ■ 

J 

) 


(mo 


1 / 

+ ■ 

— 


u. 


\ /^V0, 

V0,\ 

• (u, - Uj) 

+ 

— 

/ V "Tj 

mi y 



47r miiTij , o 

O TTT-o 


-20, + 30, + 


m,0 




Wi 


(V • Vj) — — [ —20J + 30i + 


Too 


mi9j 

m-i 


mj9. 


mo 

+u,: • 


V0i V0J 


TOi 

+ U 


V0 


1 _ ± fss, _ 2J + ”>.«j 


mo 




(V • V,) 
V0, 


m-i _ , TO,- _ , 

— (Vv, - Ui) - — (Vvj -Uj) 
9% t/j 


(— + ^) 
rrij m-j 


+ ^ _ 2Vu, • V0,) - — f ^ - 2 ) (uj • V0,) 

mo / mo \mi ' 

+u,-5[V(vi +Vj),(ui - Uj) 


(51) 
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Eq.(15) 


] Rahaman’s results 
] Huilin's results 


10 


-11 

X 10 



X 10 


Figure 2: Surfaces of scaled collision frequency in terms of two granular temperatures 9i and 
6j of our results (Eq. ID) , Rahaman’s result m and Huilin’s result na, where rrii = 
2.79e — 10 kg, rrij = 6.62e — 10 kg. Nij is scaled by d^^ninjg. 

5. Results and discussion 

We now evaluate our theory predictions by comparing with the previous 
works. In order to particularize the general results given in the the above 
120 section, we adopt the typical fluidized bed conditions used in Ref. [53]. The 
parameters include masses rrii = 2.79e — 10 kg, nrij = 6.62e — 10 kg, diameters 
di = dj = 5e — 4 m, the restitution coefficients Cij = Ci = Cj = 0.9, the granular 
temperature ranges from 1.5e — 11 kg m? js'^ to lOe — 11 kg rr? j , the volume 
fractions = 0.25, Smax = 0.638. 

125 We begin with the collision frequency iVy which presents the number of 
collisions between particles of species i and j per unit time and per unit volume. 
Nij is expected to be equal to Nji due to the symmetric character in collisions. 
So this implies that the expression of Nij can not contain the odd power of 
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Figure 3: Profiles of scaled collision frequency as a function of granular temperature Oj, where 
9i = 4e — 11 kg vn? j = 2.79e — 10 kg^ and mj = 6.62e — 10 kg. Nij is scaled by 
df-riinjg. 


B. Fig. ([^ plot iD surfaces of scaled collision frequency of Eq. (151 in terms 
130 of granular temperature of two species and comparison with previous results 
[231 [^. Nij is scaled by dfjrunjg considering that Nij vary greatly in different 
value of the Enskog factors (here, denoted by g). Three surfaces intersect at 
line 6i = 6j. Furthermore, we can observe that our surfaces are higher than 
the Huilin’s prediction and lower than Rahanman’s when 9i > 9j. This is 
135 expected as the Huilin’s profile can only be applied in energy equipartition 
systems as pointed in [26], which may underestimate the collision frequencies, 
whereas Rahaman’s hypothesis implies all the collisions happen in a plane, which 
obviously overestimates the collision frequency. When 9i < 9j , The value of Eq. 


(15) is still smaller than Rahanman’s results. Variation or error of Nij of Huilin’s 
140 results become large due to negative B. 

To have a closer examination, we plot a cross section of Fig. ([^ as a function 
of granular temperature 9j when 9i = Ae — 11 kg m? j which is illustrated in 
Fig. (§. It can be seen that our results intersect with the Rahaman’s and 
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Figure 4: Surfaces of scaled collision frequency in terms of two granular temperatures 6i and 
where (a) rrii = mj = 2.79e — 10 kg] (b) irij = 2.79e — 10 kg, m* = 6.62e — 10 kg] (c) 
rrii = 2.79e — 10 kg, mj = 6.62e — 10 kg. Nij is scaled by d^-nirijg. 


Huilin’s results at the point 9i = 9j = 4e — 11 kg rv? j which confirms that 
145 these three results are equal at 9i = 9^. Besides, we could find that Rahaman’s 
Nij decreases in the region 9i > 9j, and increases in the region 9i < 9j with 
increasing 0^, but our results show that Nij increases with increasing 9j in 
both of two regions. As the collision frequency increases with the increase of 
granular temperature, the value of Nij of the region 9j < 0^ = 4e — 11 kg rr? j 
150 is supposed to be smaller than that of point 9j = 0^ = 4e — 11 kg vr? j. So we 
can expect that our results are more reasonable than Rahaman’s. 


In Fig. Q, We plot N^ in Eq. (15) in terms of granular temperature of 
two species under various mass ratios, (a) vrii = vtij = 2.79e — 10 kg-, (b) mi = 


2.79e — 10 kg, mj = 6.62e — 10 kg-, (c) mj = 2.79e — 10 kg, mi = 6.62e — 10 kg. 
155 Surface of collision frequency Nij with mi = mj (a) is symmetric about the 
9i = 9j. If mi 7 ^ mj, the surface of collision frequency Nij tilts about 0^ = 9j 
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Figure 5: Surfaces of collision stress Pj as a function of granular temperature 6i and Oj^ 
where = 2.79e — 10 kg, rrij = 6.62e — 10 kg, and Sij = 0.9. The main figure includes three 
results: Rahaman’s, Arastoopour’s and Eq. (|28|l and the inset figure adds Huilin’s results. 


axis. We exchanged the mass rrii and the surface (b) and (c) is symmetric 


about 9i = 6j . 


Fig. (§ show the variation of the collisional stress component (zeroth-order 
160 approximation) with two granular temperatures. The main figure includes three 


results: Rahaman’s, Arastoopour’s and Eq. (28), and the inset adds Huilin’s 


results. It could be found that our results are tangent to Arastoopour’s at 
9i = 9j and are very close to Rahaman’s results. In Arastoopour’s results [2H] . 
ij reaches its maximum near the point of maximum of \9i — 9j\, i.e (1.5e — 
165 II kg rr? j, lOe — II kg rr? js^) and (lOe — II kg rr? j^ 1.5e — II kg rn? js^) 
in Fig. (j^. While our results shows the maximum of P^ij at maximum of 
both of 9i and 9j, i.e. (lOe — II kg vr? j— 11 kg rr?js^) in Fig. 

Such difference may be caused by that Arastoopour’s work is based on the 
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Figure 6: Surfaces of scaled viscosity as a function of granular temperature 6j and 6j, where 
mi = 2.79e — 10 kg, mj = 6.62e — 10 kg, Si = Sj = 0.25, Emax = 0.638, di = dj = 5e — 4 m, 
and Ei = Bij = 0.9. 


assumption that interaction between particles from two species is only at the 
170 interface, while our results are free of this assumption and hence the integration 


of Eq. (281 is complete. 

Fig. (§ plots the solids shear viscosity (Eq. ( [M| )) and its comparison with 
Huilin’s [24] and Rahman’s [25] results. Because the Enskog factor in Huilin’s 
paper is much larger than the rest {g = 58 in paper [24j . g = 1.68 in paper [25] 
175 and g ^ 5.52 in our case (Eq. §)) , we scaled Huilin’s results [23] to make these 
three viscosities can be drawn in one figure. Unlike Rahaman’s result where the 
viscosity has a minimum at 9i = 6j, our result shows that the fii is still mainly 
determined by 9i and the maximum of viscosity appears at the maximum of 9i , 
which demonstrates that the dilute viscosity increases with the increase of the 
180 granular temperature. 

The model proposed in this article allows different granular temperatures, 
mean velocities and diameters, masses, etc. In the previous discussion, the 
parameters such as the collision frequency and shear viscosity seem to be only 
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Figure 7: Surfaces of energy source 7^ ij ^ function of mean velocities Uj and Uj, where 

Si = £j = 0.25, £max = 0.638, di = dj = 5e — A m, 6i = 6j = 3.9000e — 11kg m? 
mi = 2.79e — 10 kg, mj = 6.62e — 10 kg and Sij = 0.9. 


affected by granular temperatures and masses, but not by mean velocities of 
185 particles of two species. Previous researches [211ES] are not accountable to two 
different mean velocities, either. To investigate the influence of mean velocities, 


Fig. (m) illustrates the surface of energy dissipation component yV (Eq. (50)) 


c.ij _ 

in terms of two mean velocities. It can be found that when the mean velocities 
are not equal, varies greatly, reaching its peak when is maximum and 
190 Uj is minimum. So the difference of two mean velocities should not be ignored. 
In an air-fluidized bed, the mean velocities of clusters and dispersed particles 
are usually not equal to each other. Our model is sensitive to this factor. And 
more elaborate validation needs further efforts. 


6. Conclusion 

195 In this paper, we develop a kinetic theory based model for a three-dimensional 
binary granular mixture with different masses, sizes, mean velocities, densities 
as well as granular temperatures, using standard Enskog theory. The integra- 
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tion of the inner product of relative velocity c^- and combined velocity G is 
properly treated in three-dimensional space. 

200 The computed collision frequency and particle viscosity coefficient increase 
monotonically with the increase of granular temperature of each particle species. 
Our results also show that the energy dissipation component depends heavily 
on the mean velocities of particle species, thus, difference of mean velocities, if 
any, should not be overlooked. 

205 Compared with previous work, our research enriches and consummates the 
previous theories and is more suitable for the multi-type particle theory in gas- 
solid flow systems. This work can be applied in the case that the effect of energy 
non-equipartition and unequal mean velocities is distinct, e.g., binary granular 
mixture or fluidized bed. 

210 We must remind ourselves that the granular binary hydrodynamics equa¬ 
tions, as outlined in previous sections of this article, can only deal with the 
dilute and near elastic granular systems. We argue that our model is closer to 
the real behavior of dilute granular flow in air-fluidized beds than the previous 
results. 
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